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We present results for the high-temperature series expansions of the susceptibihty and free energy 
of the g-state Potts model on a D-dimensional hypercubic lattice for arbitrary values of q. The 
series are up to order 20 for dimension D < 3, order 19 for D < 5 and up to order 17 for arbitrary 
D. Using the g — > 1 limit of these series, we estimate the percolation threshold pc and critical 
exponent 7 for bond percolation in different dimensions. We also extend the 1/D expansion of 
the critical coupling for arbitrary values of q up to order D~^. 
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I. INTRODUCTION 

Since many years the g-state Potts model l|WtJ . Il982() has served as paradigmatic system for the study of phase 
transitions and critical phenomena. Fortuin and Kasteleyn reformulated it as a random cluster model in which q 
appears as a general parameter no longer restricted to q — 2,3, . . . and where bond percolation appears as the special 
case q = I- 

Systematic series expansions have a long history l)Domb and GreenL Il974l) as a method to analyse lattice models 
of phase transitions. In the case of Potts models, high-temperature series for the susceptibility and free energy with 
q as a free para meter have been obtained for ZJ-dimensional hypercubic lattices in fMoraal, 1993) (free energy only) 
and l|Schreider and Regeii Il993) . We extend these series considerably and use them for a scan of the phase diagram 
in (g, D)-space. A detailed analysis is presented for the case of bond percolation {q = 1) on higher -dimensiona l 
hype rcubic lattices, extending the analysis of 15th order percolation series in arbitary dimensions of l|Adler et al\ . 
[199(f). 

We shortly review the model and its cluster formulation in Sect, m and the star-graph expansion technique in 
Sect, mil In the following sections we present results for general parameter q in two (Sect. Hvjl and higher (Sect. Ivjl 
dimensions. Section IVTI treats the large-dimensionality expansion of the critical coupling, and in Sect. IVlTl the case of 
bond percolation, i.e., the limit q = 1 is analyzed. We end with a short summary in Sect. IVlTll 

II. THE g-STATE POTTS MODEL 

Let G = {V, E) be an undirected graph, that is a set V of vertices together with a set E of edges defining which 
pairs of vertices are considered "nearest n eighbors" . Put a spin Si on each vertex i € V which can take q different 
values. The (/-state Potts model (see, e.g., l|Willl983) for a review) on this graph is then defined by the Hamiltonian 

H = -J 5{s,,s,), (1) 
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where S{., .) is the Kronecker delta symboL The Fortuin-Kasteleyn representation shows that the partition function 
Zg ~ } for a finite graph G is actually a polynomial in q and w — e^'^ — 1: 

(2) 

CCB 

Here the sum is over all 2l^l clusters, i.e. subsets C of the set E of edges; Nc{C) is the number of connected components 
where each isolated vertex is counted as one component and \C\ is the number of edges of C. This allows one to 
study the model for arbitrary values of q ^ {2, 3, . . .} which can no longer be represented by a local interaction of spin 
degrees of freedom of the original Hamiltonian . Of special interest in statistical physics are the limits q ^ 1 and 
q — !■ which can be shown to be related to percolation and tree percolation on the graph G, respectively. 
A two-point correlation function can be defined in the cluster representation simply by 

m,j) = ^ E 9"^=""^-'^', (3) 

where the cluster sum is restricted to those clusters Cij in which the vertices i and j belong to the same connected 
component. The susceptibility related to this correlator, 

XG{q,w) = 7^ E ^(*'-^')' W 

equals for q ~ 2,3, . . . the usual magnetic susceptibility of the g-state Potts model in the high-temperature phase (i.e., 
as long as (si) — 0). 

In the case g = 1 the weight of a cluster in the partition sum depends only on its number of edges. This case 
therefore describes bond percolation where the local bond probability is given hy p — w/{w + I). It easily follows that 
here the susceptibility x{q = ^^w) measures the average size of the percolation cluster. 

The Potts partition function Zaiq, w) is essentially equivalent to the Tutte or dichromatic polynomial of the graph 
G, see, e.g., JSokal, 200i and references therein. It encodes a lot of combinatorical properties of the graph such as the 
number of spanning trees, the number of spanning forests, the number of possible acyclic orientations, or the number 
of vertex colourings with q different colours using different colours on neighbouring vertices. The last mentioned 
number, for example, equals hm^j_^_oo ^g(<Z, e^"' — 1) = ZG{q,—l), which is the number of ground states of the 
antiferromagnetic g-state Potts model on the graph G. 

In statistical physics, on the other hand, the interest in these models is related to the existence of a phase transition 
between a disordered high-temperature (low-density in case of percolation clusters) and an ordered low-temperature 
(large w) phase if we go to infinite graphs like the Z?-dimensional hypercubic lattice G = iP . Depending on q and 
D, these transitions can be first or second order. Whereas in Z? = 2 dimensions many exact results are available 
showing that the second-order nature changes to first-order for g > 4, in I? = 3 and higher dimensions there is strong 
numerical evidence that the transition is of first order for all integer values q > 2. In the case of a second-order phase 
transition at /3c, the behavior of the susceptibility Q is characterized by the critical exponent 7: XG ~ |/3 — l3c\~'^- 



III. HIGH-TEMPERATURE SERIES EXPANSIONS 

Series expansions are an important tool for the extraction of information on statistical systems which is exact 

i up to the order of the expansion) in the thermodynamic limit. There exist different well-established methods 
Domb and Greeil Il974|) for the systematic generation of high-temperature series expansions, as for example the 
linked cluster and the star-graph method. The latter one got its name from the fact that only 1-vertex-irreducible 
graphs contribute to the series. In an attempt to study sy stems with quenched disorder via series expansions 
llHellmund and JankeL l2002albl. l2003i \2m^ iJanke et all l2005(l we used this method and developed a comprehen- 
sive software toolbox for generating and enumerating star graphs which is re-used in this work. 

The basic idea of the star-graph method is to assemble the value of some extensive thermodynamic quantity F 
on a large or even infinite graph from its values on subgraphs: Graphs constitute a partially ordered set under the 
"subgraph" relation. Therefore, for every function F{G) defined on the set of graphs exists another function Wf{G) 
such that for all graphs G 



F{G) = E ^^^9), 

g<ZG 



(5) 
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FIG. 1 Two star graphs of order 19 (left) and 17 (right) and their weak embedding numbers. Due to its topology, the result is 
complete (vahd for arbitrary dimension D) for the graph of order 17. For the graph of order 19 it is valid only up to _D = 6. 



and this function can be calculated recursively via 

Wf{G)^F{G)-Y,Wf{9). (6) 

SCG 

This gives for an infinite (e.g. hypercubic) lattice 

F{Z")^Y,{G:Z^)Wf{G), (7) 

G 

where (G : Z^) denotes the weak embedding number of the graph G in the given lattice structure l|MartinL Il974l) . 

The following observation makes this a useful method: Let G be a graph with an articulation vertex where two star 
subgraphs Gi^2 are glued together. Then Wf{G) vanishes if 

F(G) =F(Gi)+^^(G2). (8) 

An observable F for which eq. ((HJ is true on arbitrary graphs with articulation points allows a star-graph expansion. 
Then all non-star graphs have zero weight Wp in the sum eq. iQ. It is easy to see that the (properly normalized) 
free energy InZ has this property and it can be shown that the inverse susceptibility 1/x has it, too. 

The weak embedding numbers (G : Z^) of a graph G in to the hypercubic lattice are counted using a refined 
version of the backtracing algorithm by Martin ( !Martinl . ll974|) . making use of a couple of simplifications for bipartite 
hypercubic lattices Z^. Figure ^ shows two typical results. A star graph of order n can use at most [^J dimensions 
of the lattice and every embedding using m dimensions appears in Z^ with multiplicity (^) . 

For the symbolic calculations of the partition sum and susceptibility on every star graph up to a given order (ca. 
80 000 up to order 20) as polynomials in q and w, we developed a C++ template library using an expanded degree- 
sparse representation of polynomials and series in many variables. The open source library GMP is used for the 
arbitrary-precision arithmetics. Our series for the free energy In Z and the susceptibility x with g as a free parameter 
are up to order 20 for dimension D < 3, order 19 for D < 5 and up to order 17 for arbitrary D. 

In order to check our algorithms we did a careful comparison with older p ublished series data. Our se ries agree with 
the 10th order susceptibility and free energy series for arbitrary q and D of ijSchreider and Regerl [l,993) (modulo some 
misprints in their expression for In Z) as well as with the 16th order D = 2 free energy series offMoraal. 1993) but 
disagree in several places with the 11th order, arbitrary dimension free energy series in the last mentione d paper. For 
the special case of bond percolation q = 1 they agree with the 15th order arbitrary-dimension series of l|Adler et all 
^90). In the special case of the three-dimensiona l Ising model {q ~ 2, D ~ 3) our 20th order series also conform with 
the 25th order series of (|Biitera and Conil [2002) ge nerated by linked-cluster ex pansion. They also agree with the 
21th order q = 3, D = 3 partition function series of (,G!uttmann and EntineL[l994l) obtained by finite lattice methods. 
In two dimensions, much lo nger series exist for special m odels. Our 20th order susceptibility series agree with the 
q = 3 and 4, = 2 series of l|Enting and Guttmannl . l2003() . 

We close this section with some remarks about the analysis of the thus generated expansions. The estimate of 
critical parameters from a high-temperature series involves extrapolation from a finite number of exactly known 
coefficients to the asymptotic form of the function. Many such extrapolation techniques have been developed and 
tested for different series and are comprehensively reviewed in (Guttmann, 1989). These extrapolation techniques are 
not rigorous. They make some assumptions about the expected form of the singularity at the critical temperature. 
Usually, error estimates rely on the scatter of the results of extrapolations with different parameters (like [N/M] Fade 
approximants for different values of N and M). This may seriously underestimate systematic errors coming from 
wrong assumptions about the structure of the singularity. 

In order to get a reliable picture, we will take into account several criteria, such as 
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FIG. 2 Critical coupling Wc as function of g in 2 dimensions. 



• convergence of the analysis, 

• scatter of different approximants, 

• number of defective approximants, 

• agreement between different extrapolation methods. 

The basic meth ods we use are DLog-Pade approximation and inhomogeneous differential approximants (IDA) 
l)Guttmannl . ll989|) . In order to analyze confluent nonanalytic and logarithmic corrections, these methods are applied 
to suitably transformed forms of the series. T he parameters of th ese transformations are fine-tuned according to the 
criteria listed above, a technique pioneered in l|Adler et oiJ . Il99(il) . 



IV. TEST: TWO DIMENSIONS 



The Potts model in two dimensions is exhaustively reviewed in ij.Tacobsen et fl/J . l2fl0fj) . We repeat here the principal 
facts: Baxter has determined the exact free energy on 1? along the two curves 

WF = ^/q (9) 



WA = -2+^4^ (10) 

in {q, w) parameter space. The curve Q gives the locus of the ferromagnetic phase transition and is perfectly 
understood in terms of conformal field theory. Critical exponents are exactly known, e.g. 

4 + 3fc^ 

7 ^ i2fc- 6P ""^^"^^ = -2cos(7r/fc). (11) 

The curve ()lUf) is conjectured to be the locus of the antiferromagnetic transition which agrees with the known facts 
for q = 2,3,4. 

An analysis of the susceptibility series of order 20 using Fade approximants to the logarithmic derivative of the 
susceptibility (DLog-Pade approximants) is shown for the ferromagnetic case in Figs. El and |3I where around 15-25 
different [N, M] approximants are plotted for each value of q. The continuous lines present the exact results © 
and 1)11(1 . We find a quite satisfactory agreement also for noninteger values of q, however, the scatter of different 
approximants and the number of defective approximants increa ses for smaller values o f q. This may be related to the 
fact that the g — > limit is a quite intricate multicritical point ijJacobsen et a/.Ll2005(l . 
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FIG. 3 Critical exponent 7 as function of g in 2 dimensions. 



FIG. 4 Critical coupling Vc as function of q in different dimensions. 



V. HIGHER DIMENSIONS: OVERVIEW 




The explicit q- and D-dependence of our series allows us to get a fast overview over the phase structure in a large 
parameter range. Figures ^ and |5l show the results of a DLog-Pade analysis of the susceptibility series for the critical 
couphng V — {e^'^ — l)/{e^'^ — 1 + q) — w/{w + q) and exponent 7. 

Several remarks are in order: 

• At the upper critical dimension {D ~ A for q — 2, D — 6 for q — 1) the DLog-Pade analysis finds a value 
7 « 1.08 slightly larger than 1. This effect is well-known and due to the existence of logarithmic corrections. 

• For D > 3, q > 3 the phase transition is of first order and the correlation length remains finite. In analysing 
high-temperature series by Fade approximants or similar techniques, the approximant will provide an analytic 
continuation beyond the first-order transition temperature Tg into a metastable region on a pseudo-spinodal 
line with a singularity T* < Tq and effective exponents at T* . Our Vc in this region is therefore an upper bound 
of the real transition point. 



FIG. 5 Critical exponent 7 as function of q in different dimensions. 



• For D — the curve ^{q) passes 7 = 1 at g « 2.7, marking the onset of a first-order transition. This is 
consistent with shghtly smaUer estim ates in th e hterature coming from Monte Carlo computer simulations and 
field theoretic calculations (see (^GroUau et alV I2OOI ') and references therein) . It would be interesting to see if 
the crossover point can be understood as a kind of tricritical point in the (5, v) parameter space. 

Since high-temperature series alone are not sufficient to analyse first-order transitions (and much longer series exist 
for the special case of the three-dimensional Ising model q — 2, D — ^i), we concentrate our analysis in the second 
part of the paper on the q = 1 case of bond percolation. 



VI. LARGE-DIMENSIONALITY EXPANSION 

The dimension of the lattice enters the star-graph series expansion only throu gh the embedding numbe rs (eq. (0) 
and we keep it as a free parameter in our series (see Fig. 0. This allows one (jFisher and SingH ll99Cl() to extract 
a large-D expansion for the critical coupling from the susceptibility series by iteratively solving the critical point 
equation \/x{D^Vc) = in terms of Vc as series in 1/D. This gives the large-D expansion of Vc up to order _D~" if 
X is known as series in v up to order 2n. For the 0{N) mod el the first terms of the 1 /D expansion of the critical 
coupling for arbitrary N were obtained by Gerber and Fisher jGerber and Fislieil Il974j) . who also demonstrated by 
solving the large N limit that this expansion is only asymptotic. 

As far as we know, our analogous series for th e q-state Potts model wit h arbitrary q is new. For the percolation 
(g = 1) case the first 5 orders were calculated in l)Gaunt and RuskinL Il97i^ . We give the series for Vc in terms of the 
branching number a which is one less than the coordination number of the lattice and therefore for the hypercubic 
lattice a = 2D — 1, 



Vc{q,(r) 



8-3g 3(8- 3g) 3(68-31g + g2) ^ 8664 - 3798g - llg^ 
^ 2cr2 ^ 2a3 ^ 2^4 ^ 

78768 - 36714g -I- AObq^ - 50q 



12cr5 

1476192 - 685680g - 2760g2 - 551^3 



120-6 



24(7^ 



7446864 - 3524352? - 11204^2 _ 6588g3 
12^^ 



(12) 



By using cr as expansion variable we see that the first correction to the Bethe lattice result Vc = l/cr is absent 
for all values of q. Figure IHl shows plots of the critical coupling as obtained from eq. (I12f) by simply summing up all 
terms. Since these expansions are asymptotic series, for relatively small D a more careful treatment with appropriately 
adapted maximal order would actually be necessary, as discussed below for the percolation case q = 1. Of course, 
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FIG. 6 Plots of eq. 1121 (using all available terms, cf. text) for different values of q. 



in the case of first-order phase transitions, the expansion (|12|l yields the pseudo-spinodal transition point in the 
meta-stable region. 



VII. q^l: PERCOLATION 

The q ^ 1 hmit of th e Potts model describe s bond percolation with local bond probability p = w/{w + 1) = v. 
Standard scaling theory ijStaufFer and Aharonvl IT994f) describes the number of clusters of size s for large s near the 
percolation threshold Pc by 



n{s,p) s ^f{{pc-p)s''), 



(13) 



introducing the critical exponents r and cr. The q ^ I limit of the Potts model susceptibility describes the average 
cluster size and the number of clusters. The respective critical exponents are related to the scaling law eq. 113|l 

by 



7 = 
2-a 



3-T 
cr 

T- 1 



(14) 

(15) 



A throughout anaysis of percolation series in higher dimensions was done some years ago by ()Adler et al\ .\l99Cf]. 
Since we extended the series from 15 to 17 (for D > 6), 19 (for D — and 20 (in D = 3) terms, a careful re-analysis 
seems to be in place (for extens i ve per colation series expansions in two dimensions, see (jConwav and Guttmannl . ll995j) 
and, recently, ij.Tensen and Zifj . [200(^) 1. 



A. Up to the upper critical dimension 

The divergence of the susceptibility (mean finite cluster size) is expected to be of the form 

Xip) - A{p, - p)-'>[l + ai(pe -P)^^ + . . .] 

for D 7^ 6 and 

x{p)^A{p,-p)-^\\n{p,~p)\' + ... 



(16) 
(17) 



at the upper critical dimension D = 6 with 7 = 1 and S = 2/7. For D > 6, one expects a Gaussian fixed point with 
7 = 1. 



TABLE I Percolation thresholds (bond percolation on lP) and critical exponents. 
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MC = Monte Carlo, MCS — Monte Carlo, site percolation, RG = Renormalization group, HTS — High-temperature series 



We analyse the series u sing DLog-P ade approximants, inhomogeneous differential approximants and the methods 
termed Ml and M2 from (|Adler et oiJ . [T991,. 199(1'! especially tailored to take confluent nonanalytic correction terms 
into account. The method Ml uses DLog-Pade approximants to 

F{p) = {pc-p)^-ix{p), (18) 

which has a pole at pc with residue Ai — 7. For a given trial value of Pc the graphs of Ai versus input 7 are plotted 
for different Pade approximants and by adjusting pc a point of optimal convergence is searched. 

The M2 method starts with a transformation of the series in p into a series in ?/ = 1 — (1 — p/pc)^^ , and then Pade 
approximants to 

G(y) = Ai(l-y)-f lnF(2;) (19) 
ay 

are calculated which should converge to 7 as y — > 1, i.e. p pc- These methods are especially useful when taken as 
biased approximants with a given value of 7 or Ai as input. Our results are compiled in TableQ] where we quote for 
comparison also previous estimates based on Monte Carlo simulations and analyses of shorter series expansions. 

Generally the series analysis appears to be more difficult for lower dimensions. The Pc we quote for £) = 3 is an 
average over a large number (ca. 200) of Pade and IDA approximants. The general pattern in the scatter of data 
is that larger values of Pc corresponds with larger values of the critical exponent 7. The value for 7 we quote as 
our result in three dimensions is obtained by biasing the Ml and M2 analyses with the more precise pc value from 
Monte Carlo simulations. Moving on to = 4 and D = 5 dimensions, nothing spectacular happens and in particular 
we do not observe any indication of logarithmic corrections to the critical behavior in D = 5, in accordance with 
renormalizatio n group expe ctations and also with newer simulations, cf. the careful discussion of five-dimensional 
percolation in (|Fortunato et aL. 2004). 

Due to the theoretically expected multiplicative logarithmic corrections ()17|l at the upper critical dimension a special 
treatment of the case D = 6 is needed. The analysis of logarithmic corrections of the general form 

f{p)-^iPc~p)-''\Hpc~p)\' (20) 
is possible by ijAdler and Privmanlll98l[l calculating approximants for 



Fip) = (Pc - p) ln(pc - p) 



f'ip) 7 



fip) Pc-P 



(21) 



TABLE 11 Bond percolation thresholds for hypercubic lattices Z for dimensions D > 5. 
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where one expects for singularities of the form H2()|l that hnip 
find the values for pc and 5 listed in Table 



F{p) = 6. Using this method biased to 7 = 1 we 



B. The percolation threshold in large dimensions 



In llGrassbe rgeill20 03(l. Grassberg er compared his Monte Carlo results for pc with a large- 1? expansion up to order 
5 in 1/D ijGgmntgmd Ruskirl Il97q) . Since we extended this expansion up to order 9, a new comparison may be in 
order. This expansion is believed to be only asymptotic. Therefore we stop summing up terms of the expansion if the 
next term is larger than the current one. This recipe results in using only the terms up to order 7 in Z? = 5 and using 
all available terms for D > 5. We also analyzed our susceptibility series at fixed D using differential approximants 
and the Ml and M2 methods. Above the critical dimension, the Ml and M2 methods are used in a variant biased to 
7 = 1. Our results are collected in Table Hll 

For higher dime nsions our large-dimensionality expansion is in perfect agreement with Grassberger's recent Monte 
Carlo (MC) data ijGrassberg cr. 200^, whereas direct analyses of the susceptibility series in the respective dimension 
give slightly larger estimates for pc- An obvious outlier is D = 6 where the cut-off criterion for the l/£)-expansion 
apparently does not work (taking artificially one term less would significantly improve the estimate). 



C. Mean cluster number and moments of the cluster number 



Our series for the free energy In Z do not allow a reliable estimation of the critical exponent a. This phenomenon is 
already known from, e.g., Ising model series. One needs higher derivatives ^ g^^^ in order to see a singular contribution 
which shortens the series and due to the nature of the non-singular background the usual extrapolation methods give 
much worse results than, e.g., for the susceptibility series. 

On the other hand, we can determine some interesting non-universal quantities such as the mean number of clusters 
per site, 



K)= hm (N,)/\\G\l 

||G||-»oo 



(22) 



and its nth moments. Here, ||G|| is the number of sites of the graph G. Using w = qjz:^ and denoting the number of 
edges of the graph by \G\, we rewrite the partition sum as 

1 



(l-p)IGI 



^gWc(C) + |C| p\c\ 



\G\-\C\ 



implying that 



m + \c\r) 



limifg|-| Z. 



(23) 



(24) 



The mean number of active bonds per side rib — |C|/||G|| is for a I?-dimensional hypercubic lattice given by 

(rib) = pD and we get 



(ric) = lim — ^ — pD and 
9^1 oq 



||G||[((n, + nb)2)-(n, + nb)2] =lmi^(^g|^^, (25) 



TABLE III Mean cluster number and fluctuation. 



Dim. 


(nc) 


||G|| \{{n,+ntf)-{n,+nif] 


2 


0.097 9(1) 


0.161(3) 


3 


0.272 89(3) 


0.029 7(5) 


4 


0.365 494(2) 


0.007 56(2) 


c; 
o 


n 41 1 8'^9^^'^~l 


n 00"^ nA('i\ 

u.uuo uy:\^o J 


6 


0.436 327(5) 


0.00164(1) 



where / is the free energy hi Z per site. 

By using Pade extrapolation to estimate df /dq and d'^f/dq'^ at the percolation threshold Pc, we obtain the estimates 
for the mean cluster numbers and the variances eg. 1251 given in Table Hill As usual, the error estimates in Table Hill 
characterize the scatter of Pade approximants and do not include systematic errors which are presumably much larger. 

For the first number in Table IIIII the number of clusters per site (ric) in critical D = 2 bond percolat ion, 
Temperley and Lieb derived an exact expression in 1971. This expression was simplified in l)Ziff et al I Il997|) to 
(ric) = fti 0.09 8076 which affirms our analysis. For bond percolation in three dimensions, Monte Carlo simula- 

tions llLoreiiz and Z ifi. 1998 3) obt ained (tIc) — 0.272 9310(5). The fluctuations were studied in a recent Monte Carlo 
simulation l|Deng and Yanel l200(il) . Their results given for ||G|| [(n^) — (ric)^] and \\G\ \ [{ricnb) — {nc){ni,)] together 
with ||G|| [(n^) — (rif,)^] = |?5(1 — p) where z = 2D is the coordination number of the lattice and p equals Pc imply 
I |G| I [{{ric + rib)'^) — {ric + fib)^] = 0.164 45(8) for two dimensions and 0.030 71(45) for three dimensions, slightly larger 
but still consistent with our data. 

VIII. SUMMARY 

We successfully applied the method of high-temperature series expansion to g-state Potts model on hypercubic 
lattices. Modern computer facilities enable the calculation of such series while keeping some parameters (like q or 
even D) symbolic. This allows one to scan a whole parameter range. Even for the special case q — I we extended the 
known series by several terms and the results of the singularity analysis are comparable to and in good agreement with 
Monte Carlo data. By reverting the inverse susceptibility series, we derived the large- 1? expansion for the (pseudo-) 
transition point of general g-state Potts models up to order D^^. Further applications such as the study of the tree 
percolation limit q are conceivable, too. 
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